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Abstract. In this paper, we focus on Fourier analysis and holographic transforms for 
signal representation. For instance, in the case of image processing, the holographic rep- 
resentation has the property that an arbitrary portion of the transformed image enables 
reconstruction of the whole image with details missing. We focus on holographic repre- 
sentation defined through the Fourier Transforms. Thus, We firstly review some results in 
Fourier transform and Fourier series. Next, we review the Discrete Holographic Fourier 
Transform (DHFT) for image representation. Then, we describe the contributions of our 
\q \ work. We show a simple scheme for progressive transmission based on the DHFT. Next, 

we propose the Continuous Holographic Fourier Transform (CHFT) and discuss some 
■ theoretical aspects of it for ID signals. 

> 

•i-H , 

^ ' 1 Introduction 

h ; 

The techniques for signal representation composes a class of methods used to characterize the information therein 
fl2iri6l . In the area of digital image processing, a classical approach for image representation is based on unitary 
matrices that define the so called image transforms. The Discrete Fourier Transform (DFT) is the most known 
example in this class. 

Stochastic models describe an image as a member of an ensemble, which can be characterized by its mean 
and covariance functions. This allows development of algorithms that are useful for an entire class or an ensemble 
of images rather than for a single image. Fall in this category the covariance models, the 1 — D models (an image 
is considered a one-dimensional random signal that appears at the output of a raster scanner) and 2 — D models 
(Causal, Semicausal and Noncausal) lITol . 

Multiresolution representations, such Pyramid fP71 and Wavelets [20|, compose another important class in 
this field. Finally, for image understanding models, an image can be considered as a composition of several objects 
detected by segmentation and computer vision techniques fT2llT6l . 

In this paper we focus on the first category of image representation methods. We firstly review some results 
in Fourier transform and Fourier series and show some connections between them (section Next, in section |3j 
we describe the discrete Holographic Fourier Transform (DHFT) proposed in [5 1, which is based on the DFT. The 
holographic representation of an image has the property that an arbitrary portion of the holographic representation 
enables reconstruction of the whole image with details missing. This is an interesting property that motivates 
works using holographic representations for image compression [2|, multimedia systems [ 11 ], watermarking (see 



[5 1 and references therein) and representations of line images llT3l . Besides the holographic representation based 
on the DFT there are also approaches that use subsampling methods @. 

In sections |4] and |5j we show our contributions in this field. We propose the Continuous Holographic Fourier 
Transform (section HJ) and discuss some theoretical and numerical aspect of it. Then, in section |5j we develop a 
simple scheme for progressive transmission of images through the DHFT Section[6]presents our conclusions and 
some comments. The Appendix reviews interesting results in discrete Fourier analysis. An extended version of 
this work can be found in Q 

2 Fourier Analysis Revised 

In this section we review some important results about continuous Fourier Analysis. A special attention will be 
given to the connections between the Fourier Transform and Fourier Series representation of a signal. Thus, let 
us begin with usual definitions in this area. For simplicity, we restrict our discussion to the one dimensional 
case. Generalizations to higher-dimensional cases are straightforward. In this section our presentation follows the 
reference Q. 

Throughout this text, all functions / are defined on the real line, with values in the complex domain C, 
f : 5ft — ► C. Besides, they are assumed to be piecewise continuous and having the integral: 

\f(x)\ p dx, (1) 

finite. The set of such functions is denoted by L p (5ft). In this set, the p-norm and the inner product are well 
defined by: 

1 Vp 

\f(x)\ p dx\ , l<p<oo, (2) 



/oo 
f(x)g(x)dx, / l9 a 2 (»). (3) 
-oo 

Definition 1: Given a function / £ L 1 (5ft) , its Fourier Transform is defined by: 

/oo 
f(x)exp(-torjwx)dx. (4) 
-oo 

Definition 2: Given the Fourier transform /, the transformed function / can be obtained by the Inverse Fourier 
Transform, given as follows: 

/ (x) = (f- 1 /) (x) = I" 7(w) exp (2irjxu) duo. (5) 

J — oo 

Definition 3: Let / and g be functions in L 1 (JR.). Then the convolution of / and g is also an L 1 (5ft) function 
h which is defined by: 



h (x) = (f* g )(x)= f{x-y)g (y) dy. (6) 

J — oo 

Theorem 1: Any function / G L 2 (-L/2, L/2) has a Fourier series representation given by: 



+oo 



f ( x ) = a « ex p ( 2j7rx x) ' ^ 

n=— oo 

Dem: See O. 

Among the properties of the Fourier Transform, the following ones have a special place for signal processing 
techniques, according to our discussion presented in [ 1 ]. 



Property 1: If the derivative df /dx of / exists and is in L 1 then: 



dx 



(w) = /' (w) = jw/. (8) 



Property 2: f — > 0, as — > oo or — > — oo. 
For the proofs, see [7], pp. 25. 

Property 3 (Convolution Theorem): The Fourier transform of the convolution of two functions is the product 
of their Fourier transforms, that is: 

h(x) = (f*g)(x)&h{u) = f{u)g{u). (9) 

Dem: See 0. 

An important aspect for signal (and image) processing is the relationship between the coefficients a n in the se- 
ries Q an d the Fourier transform defined by the integral ©, in the case of a function / G L 2 (— L/2, L/2). In or- 
der to find this relationship, let us remember that, from the orthogonality of the functions {exp (2j7rx , n G Z} , 
it is easy to show what the coefficients a n can be obtained by the expression: 

1 [ L / 2 ( n\ 

a n = t / f( x ) ex P - 27T jx T dx. (10) 
L J-l/2 v L ' 

Now, we must observe that this expression is equivalent to: 

I f( I )= I lj(x)e, P (-2 7 r J x I )dx, (11) 

once / G L 2 (-L/2,L/2) . 
Thus, 



This result shows the connections between Fourier series and the Fourier Transform. It is used in Q to 
discussion Fourier analysis in te context of signal processing. 

3 Discrete Holographic Fourier Transform 

Optical holography technology uses interference and diffraction of light to record and reproduce 3D information 
of an optical wavefront [3 19 1. This technology has been applied for 3D displays [18| and data storage devices 
(l5i . From an arbitrary portion of an optical hologram that encodes a scene, a low quality version of the entire 
scene can be reconstructed. The quality of the recovered signal depends only on the size of the hologram portion 
used, but not of the place from where it was taken [8]. The holographic representation of images tries to mimic 
such property through image transforms. In this section, we summarize the results presented in [5 1. 
The Discrete Holographic Fourier Transform (DHFT) of an image I = I (x, y) is defined as lIBI : 



H(u,v)=F- 1 



I{x,y) e 27TjP{x > y) 



(13) 



where P (x, y) is a random phase image so that E [P (x, y) P (x, y)] = for (x, y) ^ (x,y) and P (x, y) is a 
random variable uniformly distributed over [— 1, 1] . By holographic representation, we mean that from an image 
portion H c (u, v) cropped from anywhere in the complex image H (u, v) , we can, by 2D Fourier transformation, 
get a version of / (x, y) so that the degradation is proportional to the size of H c (u, v). Clearly, if H (n, v) is 
available we get back I (x,y) as the amplitude of its 2D Fourier transform, and the same image recovery process 
will be applied over cropped parts of H (u,v). Certainly, we can take a portion of an unitary transform and 
perform the same process. However the result will be strong dependent from the place from where the portion 
was cropped. 



In order to demonstrate that the image transform defined in expression (fl3l has the holographic property, let 
us follow the development found in [5] and consider the ID discrete version of the above proposed representation 
method. The one-dimensional DHFT is: 

Af-l 

H(u)=Y l(k)e^ p ^^=e^ u \ (14) 

where {P (k)} is a set of independent and identically distributed random numbers uniformly distributed over 
[0, 1]. We will represent the process of cropping a portion of H (it) by multiplying it with a window function 
W (u), given by: 

W{u) = { ^!h°+^~J?!' (15) 
w \ 0,uf[a,a+ (L-l)\, 

where a G {0, 1, M — L} for simplicity. Now, we shall consider the question: what can be recovered from 
H c (u) = H (u) ■ W (u) by the one-dimensional Fourier transform. To answer this question we call Iyy (r) the 
obtained signal. Thus: 

M-l 

/,„ ( r ) = V H (u) W (u) -^e~ j % ur = 



V H(u) W{u)^=e 

u=0 

y I Yl( k )e^ p ^-!=e^A W(u)-j= 
M-l /M-l 
U=0 \ fc=0 



u=0 

M-l /M-l 



g J M 

M 



l ( k - r )W 



M 

M-l /M-l \ 
g E E e- j ^ u(r - k) W (u) I(k)e^ p M. 



M 

k=0 \ u=0 



Let us define g (r, k) as follows: 



Thus, by changing variable (a = u — a), we find: 



L-l 



% 1 



Thus, we finally have: 



k=0 bULC M 



j27rP(fc) 



(16) 



g (r,k) = E e~^ u(r ~ fc) W» = e -i^- fc ). (17) 

«=0 w=a 



g ( r , k) = e-^M)» ^ e -i©«(»-*). (18) 
But, from the Property A2, in the Appendix: 

V^ 1 e -j2wxu _ Lc -3(L-\)*x sincL, K x ^ (19) 

^ sincKX 
Through expressions (IT9l and (fT8l we can rewrite equation (fTBT l as: 

M-l T ■ L-K{r-k) 



'M-l 



l-im Lsinc 



L-K(r—k) 
M 



. k=0 



Msinc—^-^- 



Let us define: 



Lsmc — Vj — '- 

(r-k)-- 



(20) 



Thus, we have: 



Msinc 



M-l 



7r(r— k) 
M 



(21) 



I w ( r ) = e j ^[ a+L ^] r I I (k) e^^^L (r - k) 



(22) 

\ k=0 J 

Firstly, we should verify the result when L = M, and consequently, a = 0. In this case, expression d20l 
becomes: 



M-l 



I W (r) = e-m^> £ 1 ^ eH P ^ +h ^] H ( r - k) . 



k=0 

Besides, from the definition of <pL we can see that, when L = M, we have: 

, , \ / l,r- k = 0, 

because < r, k <M — 1. Finally, by using the result (l24l in expression d23t we find: 



(23) 



(24) 



M-l 



I W (r) = e-ifcC^)' £ / (k) eM^^] H ( r - fc) 

fc=0 

e M \ 2 rj (r) e L 



r(M-l) 
2M 



= /(r)e i2 ^ p(r) . 



(25) 



So, as expected, the whole image is recovered as the amplitude of the above result. In order to know how good 
will be the approximation of the original image generated when L < M, we can perform a statistical analysis of 
the recovered complex image. 

Following 0, we will assume that I (k) = Iq constant: 



'M-l 



/ W (r) = e -J «r— J' I Y I(k)e j27T l p ^ + ^( a+L ^)UL(r-k) 

V k=0 / 

M-l 

Jo e^ k ~ r H L (r-k). 



k=0 



k=0 



Henceforth: 



\Iw {r)\ 2 = Iw {r)Iw (0* 



'M-l 



|/ | 2 YL ^ Pa{k - T) ^L (r - k) £ e-** p *Wit> L (r - I) 



. k=0 



1=0 



M-l M-l 

\h? { E*£(r-fc) + ££ ^^Vi (r - fc) e'^^'V (r - I) 

fc^i fe,Z=0 



fc=0 

'm-i 



M-l 



\h\ 2 { E^(»-*)+EE e ^^( fe -)-^a-r)j 0L (r _ fc) 0i (r _ ^ 

fc=0 fc^Z fc,Z=0 

To answer how good an estimate is |lw (r)| for the value Iq, it is used in (5 J the following result: 
Lemma 1: Let 9q, 0\, 0m-i be independent and identically distributed random numbers uniformly chosen 
from [0,1], and the random variable: 

M-l 



V = E e j27Tdk <p(k), 



k=0 



where ip (k) is a sequence of real numbers. Thus, we have the following statistics for V: 



M-l 



E(V) = 0, E(\V 2 \)=Y<P 2 (k) 



k=0 



(26) 



o V 2 



M-l 

E E ^ 2 (*) ^ 2 (o • 

k^l k,l=0 



(27) 



M-l 



As I\y (r) is of the form e^ 27z9k (f (k), it follows from this Lemma that: 



k=0 



M-l 



i%wi 2 =i/ i 2 E^( r - fc ) 



k=0 



a 



\Iw(r)\ 2 =\I \ 2 



M-l 



\ fc^Z k,l=0 



If L = M, expression d24t implies that: 



E 

as expected from the result d25l . 
Supposing L < M, it follows: 



\Iw{r)Y 



I T |2 



1% (r) 



0, 



M-l 



I%(r)| 2 =|/ | 2 E4(r-A:) 



ifti 2 E 



M-l ■ oLnir-k) 



sine 



M 



2^{t— k) 



(28) 



k=0 smc M 



Simplifying the summation: 



M-l ■ 2Lir(r-k) 



E 

fc=0 



sine 



M 



sine 



so, returning to expression d28l we find: 



E 



If 



M 



2 -nir-k) T, ' 

M 



M-l ■ 2Lir{r-k) 



sine 



M 



k=0 



sine 



2 7r(r-fc) m 



4/n 



M 



E 



\Iw{r)[ 



L 

M 



T 2 



So the image recovered is an approximation of the original image multiplied by the factor J . 



In the previous demonstrations we recovered the original image from a window in the holographic domain 
extended by zeros. We consider now recovering the image from the window without extending it. The Fourier 
transform of the cropped sequence H(u), u £ [a, a + (L — 1)], is: 



L-l 



-j^rk 



I T (r) = FT [H («)] = £ H (k + a) -±=e 
k=o VL 

Defining the extended sequence H(u), u G [0, M — 1] and its Fourier transform as below: 

H(k) = 



(29) 



H(k + a), ke[0, L-l] 
0, ke[L, M-L], 



M-l 



Iws (v) = FT 



H(u) 



i 



lTT vu 



u=0 



e J m 



L-l 



e 3 M . 



k=0 VM 

Comparing the two expressions d29t and d30l it can be straightforward noticed that: 

(L-l)L 



(30) 



r { \ L f L \ L 2L 

I Ws(v) = S J-I T \-v , v = 0, -, - 



1 •'"> 



M 



(3D 



We must also observe that: 



a+L-l 



M0 = E h ^^m 

it— n v 



e J m 



(32) 



Hence, by changing k = u — a, we have 

L-l 



J w (r) = £ir(* + a)-L ( 



i|r(fe+a)r 



k=0 



e J m 



/L-l 

V H{k + a)-=t 



e^« flr =% (r) -e- J M ar . 



(33) 



Iws (r) = Iw (r) • e^ ar . (34) 
Substituting this result in < l3TT) . we have: 



( L \ ! M 



Considering f = jjv, follows: 



It Ijjv) = ^ T l w (v) e?M™. (35) 



Since the image is the absolute value of the holographic representation, it follows: 



(37) 

Which proves that, being \It (f)| a subsampling of |lvp (x^) |, it is a good estimate of Io( mum plying 
\It (f)| by in case of / (A;) constant. 

4 Continuous Holographic Fourier Transform 

In this section we extend some developments of section |3] for continuous signals. Therefore, we start with a 
definition for the Continuous Fourier Holographic Transform (CHFT): 

/oo 
/ (x) exp (-2irjux) exp (2vr -j'P (a;)) ds, (38) 
-oo 

where P (x) is a random phase. Following the development of section |3j we will also consider a window in 
the Fourier domain and check the quality of the recovered signal. Before this, we shall arrange the terms of the 
integral in the following form: 

/oo 
exp (-2wjcjx) [f (x) exp (2irjP (x))]dx. (39) 
-oo 

We can think about the Fourier transform of the signal h (x) = f (x) exp (2irjP (x)) , which is well defined 
if / G L 2 . Now, we use a window function (low -pass filter) given by expression: 

W( U ) = { 1 ' -*- W T*' (40) 
[ L), otherwise. 

The filtering of signal h by W can be represented, in the Fourier domain by: 

G(uj) = W (lo) H (u) , (41) 

where G (u) is the Fourier transform of the output signal and H is the Fourier transform of h, that is, the CHFT 
just defined in expression d38b . The convolution theorem applied to expression PTl renders: 



g(x) = I h (y) w(x - y) dy, (42) 

which is equivalent to: 

9i x )=l f(y)exp(2njP(y))w(x-y)dy. (43) 



A simple numerical scheme to compute expression d43t offers an interesting application of Lemma 1 . In fact, 
if we approximate that integral by the Riemann sum: 



N 



9 (x) ~ [f (y n ) w(x- y n ) Ay] exp (2vrjP (y n )) 



(44) 



n=Q 



we can , from a straightforward application of Lemma 1, estimate the / (x) through: 

/ N \ 1 / 2 

^[/(yn)^(x-y„)Ay] 2 . 

\n=0 / 

This expression is used in [ 1 1 to check the quality of the recovered signal. 



(45) 



4.1 CHFT Analysis 

The first point to be demonstrated is the fact that we can recover the input signal / (x) if we take the inverse 
Fourier transform of the CHT given b y: 

/oo 
H (w) exp (2TTjiox) du. (46) 
-oo 

If we substitute the CHT H (to) by expression (I38t we obtain: 



/oo r poo 
/ /(y)exp(-27rju;y)exp(27rjP(y)) 
-oo \_J — oo 



dy 



exp (27r jojx) dw. 



(47) 



Using development of [9| we shall introduce the kernels exp (— io 2 /n 2 ) in order to deal with convergence 
problems for this integral. Therefore, we get: 



fn(x) 



exp (27r jujx) exp 



/ (y) exp (-2irjujy) exp (2ttjP (y)) dy 
Now, we can change the order of integration and re- write this expression as: 

/oo 
/(y)exp(2vrjP(y)) 
-oo 

The inner integral has analytical solution named by k n (x — y) : 

k n (x-y) = J exp [2irjuj (x - y)] exp (~^) duJ > 



du. 



(48) 



J exp [2irjuj (x — y)] exp f ^ J 



dy. 



(49) 



(50) 



which gives: 



k n (x - y) 



cos[27rju;(x — y)] exp 



CO 



11- 



J 



sin[27rja;(x — y)] exp 



n- 



du) + 



du. 



(51) 
(52) 



From the fact that sine is an odd function and exp is an even one, it is easy to shown that the imaginary part 
will be zero and the real part of the last expression can be computed as: 



f°° ( oj 2 \ 

k n (x — y) = 2 / cos[27ru>(2; — y)l exp I I dw. 

Jo \ n 2 J 



(53) 



This integral has analytical solution that can be obtained by expression [ 10 1 : 

/ cos 6x exp (—ax 2 ) dx = — \ f— exp ( — — ^ . (54) 
Jo 2 V a V 4a / 



Henceforth, the application of this expression gives: 

[2vr (x-y)] 2 n 2 



(55) 



k n (x — y) = n-y/vrexp 
Therefore, returning to equation (l49l we obtain: 

/oo 
f(y)exp(27TjP(y))k n (x-y)dy. (56) 
-oo 

It can be proved that (see [9| for details): 

lim f n (x) = f(x) exp {2k jP (x) ) . (57) 

n— >oo 

But, we have also: 

fn (x) = H (w) exp (27r jojx) exp ( ~ ) dw. (58) 

J-oo V n ) 

It is also possible to prove that: 

/oo f u) 2 \ f°° 
H (u>) exp (2irjujx) exp I I dw = / iif (w) exp (2irjuix) dw. (59) 
.oo V n / J-oo 

From d57l and (l59b it follows that: 

/oo 
il (w) exp (2njux) dw, (60) 
-oo 

and so, as expected the whole signal is recovered as the amplitude of the inverse Fourier transform of the CHT 
5 Progressive Transmission 

In this section we propose a progressive transmission approach based on the DHFT This is performed by dividing 
the image (in the holographic domain) in several equally sized portions. Then, we take the recovered signal of a 
portion and incrementally add those ones corresponding to the other pieces. 

Formally, given W\ and Wi two window functions similar to expression (fT5b . we have: 

M-l 

I m {r)=Y J H ( u ) Wi («) -7^ e ~^ Ur > (61) 

u=0 * 
M-l 1 

I W2 (r)=J2 H («) ^ («) ^e"^ (62) 
where (u) is the DHFT defined in equation (fT3t . Thus, the linearity of the Fourier transform allows to write: 

M-l 

I W1 (r) + 4, 2 (r) = V H (u) (W 1 (u) + W 2 («)) -_ e -^ ur , (63) 

and so: 



i«i 0) + 4, 2 (r) = 1^+^ (r) 



(64) 



Observe that if W\, W2 are disjoint windows, W\ + W2 is another window function. Consequently, there is 
not changes in the spectrum of the signal inside the window. 

Thus, following expression d6lT) - d62b . each recovered portion (I m (r) and I W2 (r) above) can be considered 
as a packet to be transmitted, whose arrival order does not matter. 

6 Conclusions 

In this paper we discuss some results concerning to Fourier analysis and Holographic representation of images. 

The approach described in section [5] must be analyzed from the viewpoint of progressive transmission ap- 
proaches llT4l 151. Certainly, the need of double precision representation and the known instability of Fourier 
analysis [7| may be several limitations for such method, if compared with other ones 151 [T4l . Probably, the main 
areas of application for the Holographic transform would be in distributed words and watermarking methods. Be- 
sides a more complete theoretical analysis for CHFT must offered in order to quantify how good the estimate in 
expression d43t is for the input signal /. 
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7 Appendix 

In this section we review interesting properties in discrete Fourier analysis. 

Property Al: The circular convolution theorem states that the DFT of the circular convolution of two se- 
quences is equal to the product of their DFTs, that is, if: 

N-l 

x 2 {n) = J2h(n-k) c xi(k), < n < 1, (65) 

k=0 

where h(n — k) c is the a circular kernel (see [ 16 1) then: 

DFT [x 2 (n)] = DFT [h (n)] DFT [ Xl (n)] . (66) 
Expression d66t can be written as: 

N-1 N-1 N-l 



*2 (k) e~i> k = h{k)e-i> k ^ x x (k) , 



k=0 k=0 k=0 

If we divide this expression by we get: 

N-l , N-l „ N-l 



y X2 {k) J- e -^vk = y h{k) J_ e -3^vk y {k) e -^vk_ 
to to ViV to 

Therefore, we get the same relation for the Unitary Discrete Fourier Transform: 

UDFT [x 2 (n)] = UDFT [h (n)] DFT [xi (n)] . (67) 



Now, if we define 

xi (k) 



pngl, k = a, 
0, k ^ a, 



we have that, from expression d65b . x 2 (n) = h(n — a) c , i.e., h (n) is a circular shift of x 2 (n). Hence 

N-l 

UDFT [x 2 (n)] = UDFT [h (n)] ^ x\ (k) e~ j % vk , 

k=0 

which implies 

UDFT [x 2 (n)] = UDFT [h (n)] e~ j ^- va .png 
Property A2. Let us consider the summation: 

L-l L-l 
m=0 u=Q 

By remembering the geometric series: 

L-l 

5 = ^ a u = 1 + a + a 2 + • • • + a L ~\ 

we know that: 

S 



l-a L 



1 - a 

If we apply this result in expression dHSl we get: 



u=0 

But we know that: 



L 1 1 _ ( e -i 2 * x \ L l _ e -j2irxL e -jnxL ^ e jnxL _ e -jnxL^ 

1 — e ~j2-KX l _ e -j2nx ' e ~jirx ( e jnx _ e -j-KX\ 



E 



— j2irxi 

e J 



pjlTX _ p-j-KX 

sin (irx) = — 



sin (Lirx) 



2j 

2j ' 



Therefore, expression d69l becomes: 

L-l 



E 



— j2-wxu 

e J 



u=0 

Hence, if x = we get: 



e ^ xL 2j sin (Lttx) _ ■ (L _ 1)mB sin(L7rs) 
e~i n 2j sin (7rx) sin (to) 



L-l 



u=0 

We can also re-write expression d70l by using the sine function, which is defined by: sine (x) 
Thus, that expression becomes: 

y^ e -j2nxu = c -i(L-l)nx L7rX \ SITLC (Lirx) = Lc _ j(L _ 1)nx sijlC (LlTx) 

^— ' irx ■ sine (irx) sine (irx) 



